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ABSTRACT 

We present an analysis of quasar variability from data collected during a photometric 
monitoring of 50 objects carried out at CNPq/Laboratorio Nacional de Astrofisica, 
Brazil, between March 1993 and July 1996. A distinctive feature of this survey is its 
photometric accuracy, ~ 0.02 V mag, achieved through differential photometry with 
CCD detectors, what allows the detection of faint levels of variability. We find that 
the relative variability, 5 = cr/L, observed in the V band is anti-correlated with both 
luminosity and redshift, although we have no means of discovering the dominant rela- 
tion, given the strong coupling between luminosity and redshift for the objects in our 
sample. We introduce a model for the dependence of quasar variability on frequency 
that is consistent with multi- wavelength observations of the nuclear variability of the 
Seyfert galaxy NGC 4151. We show that correcting the observed variability for this 
effect slightly increases the significance of the trends of variability with luminosity 
and redshift. Assuming that variability depends only on the luminosity, we show that 
the corrected variability is anti-correlated with luminosity and is in good agreement 
with predictions of a simple Poissonian model. The energy derived for the hypothetical 
pulses, ~ 10 50 erg, agrees well with those obtained in other studies. We also find that 
the radio-loud objects in our sample tend to be more variable than the radio-quiet 
ones, for all luminosities and redshifts. 
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1 INTRODUCTION 

Although variability is a well known property of quasars and 
Seyfert 1 galaxies, many of its properties are still subject of 
debate. As noticed by Pica & Smith (1983), studying the 
dependence of the variability on the luminosity may allow 
to verify whether the variability is Poissonian or not, that 
is, whether the source is intrinsically multiple (like compact 
supernova remnants evolving in a nuclear starburst; e.g., 
Terlevich et al. 1992) or may be considered as a single co- 
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herent source (like an accretion disk around a massive black 
hole; e.g., Rees 1984). 

A main prediction of the sub-units model, as some- 
times the Poissonian model is called, is that the variability 
(measured in magnitudes) should be anti-correlated with 
the source luminosity. Such a behaviour was observed in 
some early studies (e.g., Uomoto, Wills & Wills 1976; Pica 
& Smith 1983) and, more recently, in the analysis of the 
samples SA 94 (Cristiani, Vio & Andreani 1990) and SA 
57 (Trevese et al. 1994), as well as in the study of the 
SGP sample (Hook et al. 1994). These three samples were 
jointly revisited by Cristiani et al. (1996), that confirmed 
the anti-correlation obtained previously. However, most of 
these studies obtained a logarithmic slope for the relation 
between variability and luminosity shallower than the value 
-0.5 that is expected in the simple Poissonian model (see 
section 4.1). The same result was achieved by Paltani & 
Courvoisier (1997), who analyzed the ultraviolet continuum 
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of all Seyfert 1 galaxies, radio-quiet quasars, and low polar- 
ization radio-loud quasars contained in the IUE database. 
However, other studies have indeed shown evidence for a 
Poissonian behaviour. For instance, Cid Fernandes, Aretx- 
aga & Terlevich (1996) analyzed the variability of the SGP 
sample with respect to the luminosity and redshift, obtain- 
ing a logarithmic slope consistent with the expected value of 
-0.5. Also, Aretxaga, Cid Fernandes & Terlevich (1997) have 
shown that the variability of the SGP sample can be well- 
reproduced by a random superposition of discrete events. 

However, even the very existence of an anti-correlation 
between variability and luminosity has been disputed. For 
instance, Lloyd (1984), Cimatti, Zamorani & Marano (1993), 
and Netzer et al. (1996) have found no dependence of vari- 
ability on absolute magnitude. Such a controversy proba- 
bly arises from the different characteristics of the samples 
studied and procedures used to analyze them. Only recently 
CCDs have been introduced in quasar monitoring cam- 
paigns, improving the accuracy of data and allowing the de- 
tection of low levels of variability (see, for instance, Borgeest 
& Schramm 1994, Netzer et al. 1996 and Giveon et al. 1999). 
Besides, most of the monitored samples present a strong cou- 
pling between luminosity and redshift, typical of flux limited 
samples. Undoubtedly, this is a major problem of the avail- 
able data sets, because it prevents one of disentangling the 
real dependence of the variability on these two parameters. 
In particular, Cristiani et al. (1996) divided their data in 
slices in luminosity and redshift, and concluded that there 
was a significant positive correlation of the variability with 
redshift. This trend had already been detected by Giallongo, 
Trevese & Vagnetti (1991), and was interpreted as resultant 
of a variability vs. wavelength dependence (see below). Cid 
Fernandes, Aretxaga & Terlevich (1996) also supported this 
hypothesis to explain the positive correlation between vari- 
ability and redshift found in their fits. 

The increase of the variability with increasing frequency 
has been observed by many authors (e.g., Cutri et al. 1985; 
Edelson, Krolik & Pike 1990; Kinney et al. 1991; Paltani & 
Courvoisier 1994; Cristiani et al. 1997; Paltani, Courvoisier 
& Walter 1998). A good example of such behaviour in a low 
luminosity AGN can be seen in Edelson et al. (1996), that 
contains multi- wavelength variability data of NGC 4151, 
monitored through the AGN Watch consortium. The pa- 
rameterization of this trend is of serious concern, because 
observations in a given photometric band sample different 
parts of the rest-frame spectrum of quasars at different red- 
shifts. Of course, such effect should be taken in to account in 
the analysis of the variability of large samples of quasars. If 
well understood, this effect can be used to constrain models 
for the nature of quasars and other forms of AGN. 

In this paper we present the first results of a mon- 
itoring of 50 quasars that has been conducted at the 
CNPq/Laboratorio Nacional de Astroffsica, in Brazil, since 
1993. The observations discussed here were done with a 0.6m 
telescope in the V band during the period between March 
1993 and July 1996. The observational strategy is based on 
differential photometry with CCD detectors. This provides 
good photometric accuracy, necessary to detect very low lev- 
els of variability. These observations will be used to address 



several issues regarding the optical variability of quasars. 
In particular, we discuss here the relation of variability with 
luminosity, including a model for the dependence of variabil- 
ity on frequency. We will also re-examine whether a simple 
Poissonian model is consistent with the observations, using 
a robust fitting technique. We will also compare the vari- 
ability properties of the radio-loud and radio-quiet objects 
present in our sample. In another paper, we will present an 
analysis of the structure function of the light-curves, useful 
to estimate the variability time-scales of quasars (Garcia et 
al., in preparation). 

This paper is organized as follows. Section 2 presents 
our sample and a description of the observations and data 
reduction, as well as the light curves of the quasars. In sec- 
tion 3 we present our estimator of quasar variability and a 
simple model for correcting the observations of the depen- 
dence of variability with frequency. In section 4 we present 
the simple Poissonian model and examine whether it is able 
to explain our observations, taking into account the correc- 
tion for the variability vs. frequency effect. In section 5 we 
compare the variability observed in the radio-loud and radio- 
quiet subsamples and, finally, we summarize our conclusions 
in section 6. 



2 THE DATA 

2.1 The Sample 

The sample of 50 quasars discussed here was drawn from 
table 1 of Veron-Cetty & Veron (1987) catalog. All quasars 
have redshifts larger than 0.15 and V-band apparent mag- 
nitudes brighter than 16.5 mag in the catalog. The former 
criterion aimed to minimize the light contribution from the 
host galaxy, whereas the latter ensured that the sample was 
bright enough to make the measurement errors small, less 
than 0.05 mag. Classifying objects with radio emission at 5 
GHz larger than 10 25 W Hz -1 as radio- loud (Kellermann et 
al. 1994), there are 35 radio-loud and 15 radio-quiet quasars 
in the sample. Among the radio-loud objects, one BL Lac 
(PKS 0537-441) was accidentally selected, because it ap- 
peared in the 1987 version of the Veron-Cetty & Veron cat- 
alog as a "normal" quasar, an error corrected in more recent 
versions of that catalog. 

The sample is presented in Table 1. The columns of 
the table give, for each object, its name, the equatorial 
coordinates a and 8 (J2000), the radio class R (L: radio- 
loud, Q: radio-quiet), the redshift z, the mean apparent and 
absolute V magnitudes (V and My), the number of ob- 
served epochs (n e ), the observer and rest-frame time cover- 
age (At bs, Atrest, in years), and the observed relative vari- 
ability (8 bs = (o"/i)obs) and the mean relative measurement 
error (8 t — e/L) in the V band. Three objects in the sam- 
ple are high polarization quasars (see table 1) according to 
Veron-Cetty & Veron (1993) or NED '; one of them is the 
blazar PKS 0537-441. The apparent V magnitudes in table 1 
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Figure 1. Distribution of the sample in (a) apparent V magnitudes, (b) redshift, and (c) absolute V magnitudes. It is also shown the 
Hubble diagram of the sample (panel 1-d). 



were computed from our photometric calibration for all but 
5 quasars (see table 1) for which we adopted magnitudes 
from Veron-Cetty & Veron (1993; see below). 

The absolute magnitudes in table 1 were computed 
assuming that the optical-UV continuum of quasars may 
be described by a power-law f v oc v a with spectral in- 
dex a—— 0.3 (Francis 1996, Peterson 1997). Besides the 
fc-correction, we have also included an empirical correc- 
tion that takes into account the presence of emission lines 
and the Lyman forest on the measured continuum (Veron- 
Cetty & Veron 1993). The reddening was corrected from 
the Ab values given by NED, assuming that Ay = 0.77 Ab 
(Allen 1976). Luminosity distances were computed assuming 
H =50 km s _1 Mpc -1 , fi = l, and A = 0. 

Figure 1 shows same features of the sample discussed 
here. Figure 1-a indicates that most of the objects in 
the sample have indeed mean apparent magnitudes around 
V = 16 (median V = 16.17). The redshift distribution is pre- 
sented in figure 1-b (median z = 0.54). Figure 1-c presents 
the Mv distribution (median My — —26.4). The peak 



ogy, under contract with the National Aeronautics and Space 
Administration. 



around Mv — —24.5 is due to radio-quiet, low redshift ob- 
jects. The Hubble diagram, i.e., Mv x z, is shown in figure 
1-d. Note how strongly luminosity and redshift are corre- 
lated in this sample. This happens because, although the 
quasars in the sample span a large range of redshifts, most 
of them have mean apparent V magnitudes around 16. 

2.2 Observations and Data Reduction 

The observations were carried out at the 0.6m Boiler & 
Chivens telescope located at the Observatorio do Pico dos 
Dias, operated by the CNPq/Laboratorio Nacional de As- 
troffsica, in Brazil, during the period from March 1993 to 
July 1996. We used basically two observational setups: one 
with the imaging camera coupled directly to a 770 x 1152 
pixels EEV CCD, and the other with a focal reducer plus a 
385 x 578 pixels EEV CCD, in order to have a field of view 
of about 7' x 11' in both cases. This field is large enough to 
allow the simultaneous observation of the quasar and some 
reference stars for the differential photometry. All observa- 
tions were made in the V band, with two combinations of fil- 
ters (GG495/2+BG18/2 or GG495+BG39/3, Bessel 1990); 
we have verified that the differential magnitudes obtained 
with these two filters are essentially the same. The quasar 



© 0000 RAS, MNRAS 000, 1-?? 



4 A. Garcia, L. Sodre Jr., F.J. Jablonski & R.J. Terlevich 



Table 1. Description of the sample and observational information. 

Quasar a (2000) <5 (2000) R< 2 ' z V 



UM 18 

PKS 0003+15 
PKS 0005-239 
PG 0043+039 
Mark 1014 
3C 57 

PKS 0232-04 
PKS 0312-77 
3C 95 
3C 94 

PKS 0355-483 
PKS 0405-12 
3C 110 
PKS 0454-22 
PKS 0537-441< 1 ) 
PKS 0637-75 
PKS 0736+Olf 1 ) 
PKS 0743-67 
PKS 0837-12 
DW 0839+18 
PG 0923+201 
PKS 0925-203 
PG 1001+05 
PKS 1004+13 
PG 1008+133 
PG 1012+00 
PG 1151+117 
3C 273 

PG 1254+047 
PKS 1302-102 
PG 1307+085 
PG 1333+176 
IE 1352+1820 
PG 1352+011 
PKS 1912-549 
PKS 2021-330 
PG 2112+059 
PKS 2115-30 
PKS 2128-12 
PKS 2135-14 
OX 169 
PKS 2145+06 
PKS 2216-03 
PG 2233+134 
PKS 2243-123< 1 > 
PKS 2251+11 
PKS 2300-683 
PG 2302+029 
4C 09.72 
PKS 2344+09 



00 05 20.2 
00 05 59.2 
00 07 56.1 

00 45 47.2 

01 59 50.1 

02 01 57.1 

02 35 07.2 

03 11 55.3 
03 51 28.5 
03 52 30.5 

03 57 22.0 

04 07 48.3 
04 17 16.7 

04 56 08.9 

05 38 50.3 

06 35 46.5 

07 39 17.9 

07 43 31.6 

08 39 50.5 

08 42 05.0 

09 25 54.6 

09 27 51.7 

10 04 20.1 
10 07 26.1 
10 11 10.8 

10 14 54.8 

11 53 49.3 

12 29 06.7 

12 56 59.9 

13 05 32.9 
13 09 47.0 
13 36 02.0 
13 54 35.6 
13 54 59.8 

19 16 39.2 

20 24 35.5 

21 14 52.6 
21 18 10.5 
21 31 35.2 
21 37 45.1 
21 43 35.5 

21 48 05.4 

22 18 52.0 
22 36 07.7 
22 46 18.2 

22 54 10.4 

23 03 43.6 
23 04 44.9 
23 11 17.7 
23 46 36.7 



+05 24 12 
+16 09 50 
-23 41 17 
+04 10 24 
+00 23 43 
-11 32 32 
-04 02 05 
-76 51 51 
-14 29 08 
-07 11 01 
-48 12 16 
-12 11 35 
-05 53 45 
-21 59 09 
-44 05 09 
-75 16 17 
+01 37 05 
-67 26 26 
-12 14 34 
+18 35 41 
+19 54 05 
-20 34 51 
+05 13 01 
+12 48 57 
+13 04 13 
+00 33 37 
+11 28 31 
+02 03 09 
+04 27 35 
-10 33 20 
+08 19 51 
+17 25 14 
+18 05 19 
+00 52 50 
-54 54 47 
-32 53 36 
+06 07 44 
-30 19 10 
-12 07 04 
-14 32 55 
+17 43 49 
+06 57 40 
-03 35 36 
+13 43 56 
-12 06 51 
+11 36 40 
-68 07 37 
+03 11 47 
+10 08 18 
+09 30 46 



L 
L 
L 

Q 
Q 

L 
L 
L 
L 
L 
L 
L 
L 
L 
L 
L 
L 
L 
L 
L 

Q 

L 



L 

Q 

L 



Q 

L 
L 

Q 

L 
L 
L 
L 
L 
L 

Q 

L 
L 
L 

Q 

L 
L 



1.899 
0.450 
1.410 
0.385 
0.163 
0.669 
1.438 
0.223 
0.614 
0.962 
1.016 
0.574 
0.773 
0.534 
0.896 
0.654 
0.191 
1.510 
0.200 
1.272 
0.190 
0.348 
0.161 
0.240 
1.287 
0.185 
0.176 
0.158 
1.018 
0.286 
0.155 
0.554 
0.977 
1.121 
0.398 
1.465 
0.457 
0.980 
0.501 
0.200 
0.213 
0.999 
0.901 
0.325 
0.630 
0.323 
0.512 
1.052 
0.432 
0.673 



16.29 
15.62 
17.32 

15.89( 3 ) 

15.69( 3 ) 
16.14 
16.26 

16.10 (3) 
16.11 
16.45 

16.38< 3 ) 
14.78 
16.20 
15.80 

16.48< 3 ) 
16.09 
16.31 
16.34 
16.52 
16.54 
15.58 
16.65 
16.41 
14.98 
16.27 
15.90 
16.10 
12.76 
16.24 
14.82 
15.12 
16.06 
16.22 
18.51 
16.76 
16.30 
15.84 
16.61 
15.40 
15.83 
16.04 
16.07 
16.52 
16.50 
16.61 
15.85 
17.40 
15.95 
15.82 
16.28 



M v 

-28.64 
-26.42 
-27.08 
-25.80 
-24.07 
-26.70 
-28.16 
-24.48 
-26.56 
-27.09 
-27.24 
-27.79 
-26.90 
-26.58 
-26.10 
-26.84 
-23.86 
-28.40 
-23.73 
-27.62 
-24.55 
-24.82 
-23.32 
-25.64 
-27.95 
-24.15 
-23.83 
-26.91 
-27.38 
-26.21 
-24.54 
-26.41 
-27.30 
-25.36 
-24.97 
-28.20 
-26.26 
-26.99 
-26.87 
-24.39 
-24.43 
-27.53 
-26.86 
-24.84 
-26.18 
-25.48 
-24.86 
-27.84 
-26.12 
-26.64 



7 
6 

5 
9 
6 
6 

7 
8 
7 
7 
8 
7 
5 
5 
8 
9 
8 
8 
9 
9 

10 
10 
10 
10 
7 
8 
9 

10 
10 
11 
6 
10 
5 
8 

17 
16 
11 
10 
10 
9 

10 
9 
7 
11 
9 
8 

10 
8 
6 
7 



At obs 

(yr) 

2.74 
2.74 
2.74 
2.82 
2.02 
1.95 
1.94 
3.40 
1.94 
1.94 
3.25 
1.94 
1.94 
2.18 
3.32 
3.47 
2.41 
3.47 
2.60 
2.41 
2.60 
2.60 
2.59 
3.14 
3.14 
3.14 
3.14 
3.14 
3.14 
3.22 
0.41 
2.09 
0.33 
3.22 
3.16 
3.82 
3.22 
2.82 
2.82 
2.09 
2.67 
2.67 
2.09 
3.22 
2.82 
2.67 
3.82 
2.75 
2.75 
2.75 



^■t-rest 

(yr) 

0.95 
1.89 
1.138 
2.03 
1.74 
1.17 
0.80 
2.78 
1.20 
0.99 
1.61 
1.23 
1.10 
1.42 
1.75 
2.10 
2.02 
1.38 
2.16 
1.06 
2.18 
1.92 
2.23 
2.53 
1.37 
2.65 
2.67 
2.71 
1.56 
2.50 
0.35 
1.34 
0.17 
1.52 
2.26 
1.55 
2.21 
1.42 
1.88 
1.74 
2.20 
1.34 
1.10 
2.43 
1.73 
2.02 
2.53 
1.34 
1.92 
1.64 



(V band) 
0.038 
0.103 
0.017 
0.056 
0.084 
0.060 
0.060 
0.122 
0.258 
0.023 
0.027 
0.074 
0.028 
0.106 
1.483 
0.188 
0.547 
0.124 
0.212 
0.153 
0.093 
0.100 
0.066 
0.144 
0.027 
0.074 
0.084 
0.121 
0.057 
0.089 
0.063 
0.029 
0.011 
0.066 
0.187 
0.066 
0.035 
0.089 
0.068 
0.237 
0.140 
0.144 
0.049 
0.041 
0.106 
0.050 
0.120 
0.053 
0.133 
0.056 



(V band) 
0.022 
0.016 
0.018 
0.047 
0.048 
0.017 
0.019 
0.049 
0.016 
0.015 
0.047 
0.015 
0.016 
0.017 
0.124 
0.014 
0.033 
0.030 
0.012 
0.014 
0.012 
0.018 
0.015 
0.012 
0.014 
0.015 
0.012 
0.016 
0.010 
0.217 
0.173 
0.117 
0.043 
0.059 
0.022 
0.013 
0.016 
0.017 
0.031 
0.040 
0.032 
0.023 
0.068 
0.056 
0.084 
0.090 
0.046 
0.011 
0.079 
0.014 



Notes: (1) high polarization quasars; (2) radio class: radio-loud (L) or radio-quiet (Q); (3) objects for which the V magnitudes were 
taken from Veron-Cetty & Veron (1993). 



fields were observed with two successive exposures, each of 
them of 5 minutes, to minimize cosmic ray contamination. 
This integration time allow us to achieve a photometric ac- 
curacy of a few hundredths of magnitude for most of the 
objects in the sample. 



The data were reduced using IRAF I standard pack- 



II IRAF is distributed by National Optical Astronomy Observa- 
tories, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under contract with the National 
Science Foundation. 
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ages for CCD reduction and photometry. The images were 
first bias and flat-field corrected using the package ccdred. 
The magnitudes were measured from aperture photometry 
through the phot task in the daophot package. In each epoch 
(corresponding to each observational run) we extracted mag- 
nitudes within several aperture radii, and then chose the 
radius which minimized the photometric errors. 

Differential photometry is done with respect to refer- 
ence stars observed in the same field as the quasars. The 
procedure we devised to select these stars for each field is 
as follows. First, we selected the brightest isolated stellar 
sources as possible reference stars for that field. After, we 
computed magnitude differences between pairs of stars, se- 
lecting as reference the star that presented the lowest mag- 
nitude dispersion relative to the others during all epochs. 

By analyzing multiple observations of several fields, we 
detected a median underestimation of the photometric errors 
given by phot by a factor of 1.73, similar to that obtained 
by Gopal-Krishna, Sagar & Wiita (1995), of 1.75. After in- 
tensive study of the possible causes of this underestimation, 
we concluded that it could be originated by random fluc- 
tuations of the PSF over the image, probably induced by 
flat field errors. In fact, large scale fluctuations produced by 
the dome flat fields used in the reduction may well increase 
the errors to the observed values. To overcome this problem, 
we decided to adopt a minimum error of 0.01 mag for each 
individual magnitude extraction, since this is the expected 
photometric error of a V =16 object with the observational 
setup used in the observations. 

We performed magnitude calibrations by observing 
standard stars. These observations were done in 7 nights, 
during runs in March, May, July and August of 1996, and 
May 1997. The standard stars were selected from Graham 
(1982) list, and observed for different air masses every two 
hours during the calibration nights. As the observations were 
conducted at only one band (V), we have taken into account 
only the extinction coefficient of first order, associated with 
the airmass, to determine the apparent magnitudes of the 
reference star in each quasar field. These observations al- 
lowed us to calibrate the magnitudes of 45 quasars. The re- 
maining 5 quasars not observed during the calibration nights 
(see Table 1) were assumed to have average apparent mag- 
nitudes equal to those reported in the Veron-Cetty & Veron 
(1993) catalog. The photometric errors achieved in our sur- 
vey are relatively low (the median value of 5 e is 0.019), what 
allows the detection of low levels of variability. 

All the observations made within 5 days were grouped 
in a single observation- an epoch. The entry n e in Table 1 
refers to these grouped observations. The median number of 
epochs per quasar is 8, and the median time interval covered 
by the observations is 2.7 and 1.7 years, in the observer and 
quasar rest-frame, respectively. Table 2 summarizes some 
features of the sample, presenting the median values of most 
of the entries in Table 1, for all objects of the sample as well 
as for the radio-loud and radio-quiet sub-samples. 



Table 2. Median values of some quantities of the sample 
all sample radio-loud radio-quiet 



number 


50 


35 


15 


z 


0.544 


.614 


.385 


V 


16.17 


16.26 


16.06 


M v 


-26.42 


-26.64 


-25.36 


n e 


8 


8 


9 


At obs (yr) 


2.74 


2.74 


2.82 


At re3 t (yr) 


1.74 


1.73 


1.74 


Sobs 


.079 


.106 


.057 


Se 


.019 


.019 


.016 


.3 Light 


Curves 







The differential light curves of the quasars in our sample 
are shown in figures 2 to 8. Each data point (AV = quasar 
magnitude minus reference star magnitude) corresponds to 
the average of all the observations collected in each epoch, 
and the error bars are the corresponding photometric errors. 
The observed light curves present a rich variety of forms. 
Some objects show a monotonical increase or decrease in lu- 
minosity (e.g. 3C 273, PKS 2128-12, 4C 09.72), while others 
seem to exhibit a pulse-like pattern (PKS 1912-549). How- 
ever, the time-scale covered by the observations (and also the 
sampling rate) is not large enough to allow an appropriate 
discussion of the form of the light curves. Since we continue 
monitoring this sample, we postpone this discussion for a 
future paper, where we will discuss the data collected over 
a larger time interval. 



3 MEASUREMENTS OF VARIABILITY 

3.1 Variability Index 

The variability analysis in this paper will be done in terms of 
luminosity, instead of magnitudes. This procedure makes the 
comparison with model predictions more direct. We estimate 
the variability of a quasar with a "robust" index similar to 
that used by Hook et al. (1994): 



1 



O~obs = 



1I> 



L\ 



(1) 



where n e is the number of epochs, Li is the quasar lumi- 
nosity in the V band at the i— t h ep och, and L is the mean 
quasar luminosity. The factor y7r/2 assures that for a Gaus- 
sian distribution of observations {Li}, this index is also an 
estimator of the standard deviation of the distribution. Sev- 
eral other estimators of variability could be used (e.g., the 
ordinary standard deviation or the light-curve amplitude), 
but the results described below are, to a large extent, inde- 
pendent of the variability index employed in the analysis. 

The observed variability depends not only of the quasar 
variability, but also of the observational errors, and the esti- 
mation of the intrinsic variability a requires some hypothesis 
on the distribution of {Li} and their observational errors. 
If we assume them both Gaussian distributed, the observed 
and intrinsic variabilities are related by 



2 
°~obs 



2 . 2 
= a + £ 



(2) 
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Figure 2. Differential light curves of the sample. AV is the difference between the magnitudes of the quasar and the corresponding 
reference star. 



where e is the mean square photometric error in luminosi- 
ties. 

Sometimes it is more useful to work with relative vari- 
ability, S = a j L, than with a or a t s directly. Note that S is 
directly related to the variability measured in magnitudes, 
Om — 0.921 a/L. The relative photometric error is 8 € — e/L. 



We will hereafter designate the observed relative variability 
by 5 bs ■ In figure 9 we display S oos as a function of absolute 
magnitude and redshift. It is also shown the median variabil- 
ity values and the corresponding quartiles within bins with 
equal number of quasars (10). The two objects with largest 



© 0000 RAS, MNRAS 000, 1-?? 



Optical Monitoring of Quasars 7 



3C 95 



< 1.8 L 



1.6 



3.6 




3C 


94 






; ' 


1 


1 


' ' 


1 : 


3.5 

> 

o 3.4 


*I 




i 


i 


1 1 1 1 1 1 1 1 1 1 1 

m 
m 
m 


3.3 


- 


, 


1 


i i 


-_ 







9400 


9600 9800 


10 4 






9400 


9600 9800 


10 4 








T (JD2440000 + ) 






T (JD2440000+) 




2.9 






PKS 0355-483 










PKS 0405-12 






1 




1 1 z 




-0.1 


-' 


1 1 | 1 1 1 | 1 




2.85 


- 












- 




1 1 1 1 1 1 1 1 II 1 1 

m 
m 

m 


> 2 - B 
^2.75 

2.7 


h^ 




* i^ 


iZ 




-0.2 
> 

< 

-0.3 


Li* 


i 


2.65 




i 








-0.4 


■ i i i i 


iii 






9500 


10 4 


1.05x 


10 4 


9400 


9600 9800 


10 4 








T (JD2440000 + ) 






T (JD2440000+) 










3C 110 






2.8 




PKS 0454-22 




3.8 






1 ' 


1 - 










3.7 
> 
< 

3.6 


- i 
: i 






V 




2.7 
> 2.6 

<3 

2.5 


S 

: i 




i,iiii|iiii|ii 

H 


3.5 


— i 




1 


1 - 




2.4 


i i i 1 i 


1 1 1 1 1 1 1 1 1 


i 1 i i 



9400 9600 9800 

T (JD2440000 + ) 

PKS 0537-441 



10 4 



3 r— r- 

2 - 

> 

< i _ • 

- 



9400 9600 9800 

T (JD2440000+) 

PKS 0637-75 



T 



I I I I I I I 

9500 10 4 

T (JD2440000 + ) 




1.05xl0 4 



9500 10 4 

T (JD2440000+) 



1.05xl0 4 



Figure 3. Differential light curves of the sample - cont. 



variability in figure 9 are the high polarization quasars PKS 
0537-441 and PKS 0736+01. 

The data shown in figure 9 are consistent with a de- 
crease of the observed variability with both luminosity and 
redshift. In fact, the Spearman rank-order correlation co- 
efficient is 0.37 and -0.35, for the S t, s x My and 8 t s x z, 



respectively. The corresponding probabilities for the null hy- 
pothesis of uncorrelated data sets is small, 7.7 x 10~ 3 and 
1.4 x 10~ 2 , respectively, indicating that these trends are sta- 
tistically significant. Note that the strong correlation be- 
tween My and z present in our sample (figure 1-d) preclude 
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Figure 4. Differential light curves of the sample - cont. 



us of knowing what is the real dependence of the variability 
with respect to both My and z. 

The quasars in our sample span a large range in red- 
shift; consequently our observations sample different parts 
of their rest-frame spectra. This makes a direct comparison 
of the variability of objects somewhat meaningless, unless a 



correction for the increase of variability with frequency dis- 
cussed in the introduction is applied. Next section presents 
a model for such a correction. 
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Figure 5. Differential light curves of the sample - cont. 



3.2 Variability x Frequency 

Let us suppose that the spectral energy distribution of a 
quasar in the optical-UV region of the spectrum is a power- 
law, /„ = /*(^/i/») c ", where v is measured in the quasar 
rest-frame. The increase of the variability with frequency 
can be interpreted as a hardening of the spectrum associated 



to an increase of luminosity. Following Giallongo, Trevese & 
Vagnetti (1991), we assume that the spectral index increases 
when the source brightens, keeping the flux unchanged in 
some near-IR rest-frame frequency v* (Cutri et al. 1985). 
Edelson, Krolik & Pike (1990) have found strong correlations 
between auv and the UV flux for the CfA Seyfert 1 galaxies. 
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Figure 6. Differential light curves of the sample - cont. 



The same kind of correlation is exhibited by the radio-loud 
quasar 3C273 (Paltani & Curvoisier 1994). We model these 
spectral changes as 



a = a + 7o log I =^- 



(3) 



where a is the mean spectral index and vo denotes a ref- 
erence frequency (e.g., a filter) in the observer rest- frame. 
Despite some observational evidence that a may be corre- 
lated with the quasar average luminosity, we will assume 
that it is constant and equal to -0.3, a value appropriate for 
the optical and UV parts of the spectrum (Francis 1996, Pe- 
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Figure 7. Differential light curves of the sample - cont. 



terson 1997). This value is also consistent with that we have 
adopted for the fe-correction (section 2.1). In Eq. (3), F VQ is 
the average continuum flux integrated in a band centred on 
vq, i.e., F VQ ~ /„Ai/, where Av is the filter bandwidth. 70 
is a parameter that depends only of vq. Note that the rest- 
frame frequency that is actually measured by a filter with 



effective frequency v$ is v = (1 + z)v§, where 2 is the quasar 
redshift. 



From Eq. (3), a small change 8F VQ in the quasar con- 
tinuum flux produces a change in spectral index given by 
8a — 7o(loge)<51nK, . Since 
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Figure 8. Differential light curves of the sample - cont. 
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Figure 9. Observed dependences of the variability on luminosity 
and redshift for our sample. Crosses represent the median val- 
ues over bins of equal number of objects (10). The vertical bars 
correspond to the quartiles of the observed variability inside the 
bins, and the horizontal bars give the range in luminosity of the 
objects in each bin. 



Shift, = In 



(B 



Sa, 



(4) 



we have that the monochromatic continuum flux of a quasar 
varies with frequency as 



SJv 



=- = 70 log — -=-* 
f„ \"*/ F vn 



(5) 



Let us now consider the relative variability of the flux F V1 , 
measured within a band centred on a certain frequency i>\, 
as 

5F V1 _ 5f v 

F V1 ~ /„ 



S V1 



(6) 



where, now, v — (1 + z)v\ is the rest-frame frequency. Note 
that in the above formula we did not take into account the 
line-emission correction used before to compute the absolute 
magnitudes of the quasars. In fact, for objects with redshifts 
lower than 2.0 (our case) this correction implies differences 
of at most ~12% in Sv (i.e., 8 measured in the V band), so 
we decided to neglect such effect. Then, it is easy to verify 
that the relative variability of a quasar measured in two 
filters with effective frequencies V\ and v? are related as 



6vi _ lQg((l +z)Vi/u t ) 

5v2 log{(l + z)wt/v*) 



(7) 



We have estimated A* = c/v, « 11400 A from data 
of Edelson, Krolik & Pike (1990), by fitting Eq. (7) for 
Ai = 1450A and A2=2885A; we have used all objects stud- 
ied by those authors with the exception of NGC 7469, that 
seems to present a different variability pattern. Using data 
for 3C273 from Paltani & Curvoisier (1994), we obtain a 
similar value, A* ~ 11250A. 

This very simple model for the frequency dependence 
of variability is indeed able to explain the variability be- 
haviour of the nucleus of the Seyfert galaxy NGC 4151. For 
this exercise, we have used data collected in several spectral 
bands by Edelson et al. (1996). Firstly, we removed from 
the observed fluxes the light contribution of the host galaxy, 
accordingly to the recipe given in the paper. After, we fit- 
ted Eq. (7) to the relative variability as a function of the 
frequency, obtaining A* ~ 8400A for data in the optical 
and UV. The fit of the model to the data is very good, as 
shown in figure 10, and has a linear correlation coefficient of 
-0.92. Surprisingly, including the data point corresponding 
to the intermediate X-ray emission (1.5 keV), the correla- 
tion becomes even stronger (r=-0.99) and A* increases to ~ 
9500A. This last result, however, may be somewhat fortu- 
itous because we did not take into account any correction for 
the X-ray absorption. Anyway, our variability model seems 
to fit very well the overall variability pattern observed in 
NGC 4151. In what follows, we adopt A* = 1 /im as the 
pivot wavelength, since the results are not strongly depen- 
dent on the value of A* (see section 4.2). Note that near 1 
/jm there is also a minimum in the spectral energy distribu- 
tion of quasars, separating the bump that extends towards 
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A (A) 



Figure 10. Relative flux variation (in percentage) of the NGC 
4151 nucleus (from data of Edelson et al. 1996) as a function of 
the logarithm of wavelength. The straight line represents a simple 
(unweighted) linear fit. 



the infrared (attributed to thermal emission by dust) from 
the optical-UV Big Blue Bump. 

Let Sobs be the relative variability observed in the V 
band and Srest the relative variability that we would observe 
in the V band if the quasar was at z — 0. Since the V band 
in the observer frame corresponds to a frequency (1 + z)vy 
in the rest-frame of a quasar at redshift z (where vy is the 
V band characteristic frequency), we have that 



log(fv7^* 



Sobs \og((l + z)i>v/i>*) 



(«) 



Note that the relative variability is independent of the k- 
correction. 

These expressions can be used to transform the ob- 
served variability of the quasars of our sample to rest-frame 
variability. For instance, Eq.(8) can be rewritten numerically 
as 



J rest 



0.260 



log (1.818(1 + «)) 



(Sob 



(9) 



where we have included the effect of observational errors in 
quadrature, accordingly to Eq. (2). 

The effect of the correction is to reduce the observed 
variability. For the most distant object in our sample 
(0=1.899) this reduction is 36%. 

The application of this correction decreases slightly the 
scattering in the relations between relative variability and 
absolute magnitude or redshift. The Spearman correlation 
coefficients (rs) calculated with data corrected by the mea- 
surement errors and frequency dependence are 0.38, for 
8 2 e3t x My, and -0.36, for 8 2 e3t x z, to be compared with 
the values before the frequency dependence correction, 0.25 
and -0.22, respectively. 



between their variability and luminosity can be viewed as 
an evidence favouring the sub-units model. In the Starburst 
model (Terlevich et al. 1992), for example, the variability 
is produced by supernova explosions and the interaction of 
their ejecta with a high-density circumstellar medium. The 
superposition of random individual pulses would be respon- 
sible by the Poissonian character of this model. Even not 
being traditionally faced this way, the accretion disc model 
itself could accommodate a Poissonian behaviour, as long as 
the instabilities that would cause the variability act as inde- 
pendent events (Cid Fernandes, Aretxaga & Terlevich 1996). 
The same may be expected in models where the variability 
is produced by external effects, if these too are independent 
of each other. 

Let us consider here the Poissonian model in its sim- 
plest version (e.g., Paltani & Courvoisier 1997 or Cid Fer- 
nandes 1995). The light curve of a quasar is described as 
a superposition of random pulses over a component of con- 
stant luminosity L c , and depends only on three parameters 
related to the pulses: their rate (n), energy (E), and effec- 
tive time-scale (t p ). If l(t) is the typical luminosity profile of 
a pulse, we may define 



E 



and 



l(t)dt 



[Jl(t)dt] 



(10) 



(11) 



Jl(t) 2 dt 

Then, it can be shown that the average luminosity of a 
quasar is 

L = L c + L var = L c + nE (12) 

and that its variance is given by 

2 En EL var 



' a 



' ii 



(13) 



where L var is the variable component of the luminosity. 

Let f V ar be the fraction of the average luminosity due to 
the variable component: f var = L var /L. Then, the variance 
of the quasar luminosity can be written as 



E n Ef var L 



to 



' D 



and the relative variability is 

X _ a _ I fvarE 

~ L ~ \ L U 



(14) 



(15) 



Therefore, if fvar, E and t p are independent of the lumi- 
nosity, a and 5 must be proportional to L 1 ' 2 and L^ 1 ' 2 , 
respectively. We will call this scenario as the simple Pois- 
sonian model. In this model, quasars differ from each other 
only by their rate of events, n. 



4 IS THE VARIABILITY POISSONIAN? 

4.1 A Simple Poissonian Model 

The non-periodicity of pulses in the light curves of quasars 
and Seyfert galaxies, as well as the observed anti-correlation 



4.2 Testing the Poissonian Hypothesis 

As mentioned before, due to the strong correlation present 
in the luminosity vs. redshift diagram, we are unable to ob- 
tain with our sample the joint dependence of the variability 
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on luminosity and redshift. So, in order to test the simple 
Poissonian model, we will neglect any dependence of vari- 
ability on redshift. This assumption is justified by results of 
some previous variability studies, that have not found any 
significant dependence of the variability vs. luminosity anti- 
correlation with redshift (e.g., Hook et al. f994; Cristiani et 
al. 1996). 

In order to verify whether our data is consistent with 
the Poissonian scenario, we model the corrected relative vari- 
ability Srest as 

Srest = -jj (16) 

where a and b are free parameters. In the simple Poisso- 
nian model the expected value of 6 is 1. We will work with 
8r e3t instead of S res t because for some objects the observed 
variability is smaller than the observational errors, leading 
to negative values in the right side of Eq.(9). Our approach 
avoids discarding these data points. 

The parameters a and b were obtained with a "robust" 
version of \ 2 - Let us define the residual m of the variabil- 
ity of the i-th quasar, weighted by its relative observational 
error, as 



{51 



a/L\ 



(17) 



In traditional x 2 , the parameters a and 6 may be determined 
by minimizing the mean value of the set {« 2 }. Here we esti- 
mate these parameters by minimizing the median of {|«j|}. 
This approach is less sensitive to the variability behaviour 
of objects with extreme properties, like the high polariza- 
tion quasars and the objects with very low variability. The 
uncertainties in the parameters were obtained by bootstrap 
re-sampling of the observed data sets (Feigelson & Babu 
1992). 

The results of the application of this method to our 
data are shown in Table 3. We have obtained a — 0.05 ± 
0.25 and b= 1.20 ± 0.58 (for L v in units of 10 10 L Q ). The 
large values of the errors reflect the large spread of quasar 
variability at a given luminosity, as can be seen in figure 
9 or in figure 11, where we plot Srest as a function of the 
luminosity of the objects, both in linear and logarithmic 
scales. The continuous line represents the best fit for all 
the sample. From the bottom plot in figure 11, it can be 
seen that excluding from the analysis objects for which <5 2 est 
is negative would lead to a shallower slope (b =0.69±0.36 
instead of b =1.20±0.58) for the variability versus luminosity 
relationship. 

A very interesting result is that the best-fit value of b 
is always very near the value b — 1 expected in the simple 
Poissonian model. As shown in Table 3, even for the radio- 
quiet and radio-loud sub-samples we have obtained values of 
b consistent, within the errors, with the simple Poissonian 
model. The dashed line in figure 11 represents the best-fit 
solution for the simple Poissonian model. It is worth not- 
ing that the a and b values obtained from this analysis are 
rather insensitive to changes in the pivot wavelength A* . For 
instance, modifying A* by 1000A leads to changes much less 
than 1<t in the best-fit values of a and b. 

If we assume that the quasar variability may indeed 



be described by the simple Poissonian model, then the only 
free parameter now, a, can be estimated by minimizing the 
median of {|ui|} with 6=1. These results are also shown 
in Table 3. In this case the parameter a can be used to 
constrain the energy of the individual pulses, E, because 



E 



a x t p 

Jvar 



(18) 



Assuming that f var is in the range 0.5-1.0 (Cid Fernandes, 
Aretxaga & Terlevich 1996 adopted a lower limit of 0.3), 
and a typical pulse time-scale t p of 1.5 to 3.0 years (see, 
for example, the structure function analysis of Hook et al. 
1994, and the discussion in Cid Fernandes, Aretxaga & Ter- 
levich 1996), we get individual energy pulses falling in to 
the interval 2.5 x 10 49 — 1.5 x 10 50 erg in the rest-frame V 
band. We can compare these results with those obtained in 
the analysis of Cid Fernandes, Aretxaga & Terlevich (1996). 
For quasars, we may assume that Ly ~ 0.66 x Lb (for 
a = —0.3). Then, the energy range found by those authors 
is 1.0 x 10 50 <, E <, 4.1 x 10 51 erg in the V band. This 
energy interval is not inconsistent with our results. A recent 
study of the bright quasar 3C273 by Paltani, Courvoisier & 
Walter (1998) indicates that its pulse energy in the optical- 
ultraviolet spectral range is ~ 3 x 10 erg. The correspond- 
ing energy in the V band is about one order of magnitude 
smaller, that is also consistent with our results. 
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Table 3. Parameters obtained from the robust fits, for all data 
and for the radio class sub-samples. We present results assuming 
6=1 (simple Poissonian model), as well as for the case where a 
and b are free parameters. The symbol # refers to the number of 
data points used in the fit. The parameters were calculated with 
V luminosities in units of 10 10 Lq. 



All Sample 


a 
0.050±0.254 
0.030±0.017 


b 

1.20±0.58 

1 


number 
50 

50 


Radio-Loud 


0.130±0.395 
0.040±0.054 


1.43±0.50 
1 


35 
35 


Radio-Quiet 


0.020±0.150 
0.014±0.010 


1.22±0.81 
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Figure 11. Rest-frame variability as a function of luminosity 
for our sample. The continuous line represents the best-fit result 
assuming a and b as free parameters, whereas the dashed line is 
the best-fit for the simple Poissonian model (b = 1). The fittings 
correspond to the "all sample" parameters in tabic 3. 



5 VARIABILITY OF RADIO-LOUD AND 
RADIO-QUIET QUASARS 

Radio-loud and radio-quiet quasars (hereafter in this sec- 
tion RL and RQ, respectively) present almost the same con- 
tinuum spectral characteristics, with the exception of their 
radio emission, where they differ by about 3 orders of mag- 
nitude. It has been observed a trend for RL be more vari- 
able than RQ (Pica & Smith 1983), what can be due to 
a contribution of the collimated components associated to 
the non-thermal emission (Netzer et al. 1996). Some studies 
have suggested the existence of an intermediate class, where 
the small scale jets observed in some RQ would be directed 



Figure 12. Median rest-frame variability as a function of ab- 
solute magnitude and redshift for radio-loud (filled circles) and 
radio-quiet (open circles) quasars. The horizontal lines give the 
interval of My or z corresponding to each bin. 



along the line-of-sight (Miller, Rawlings & Saunders 1993; 
Kellermann et al. 1994; Blundell & Beasley 1998). 

The median values for the variability of RL and RQ, cor- 
rected for the variability-frequency effect are, respectively, 
Srest equal to 0.063 and 0.019. Consequently, for our sam- 
ple, the rest-frame relative variability of RL is about 3 times 
larger than for RQ. About the same factor is found for the 
parameter a in the simple Poissonian model, as shown in 
Table 3. However, RL in our sample tend to be more lumi- 
nous and have higher redshift than RQ: the median My is 
-26.6 and -25.4 for RL and RQ, respectively, and the cor- 
responding median redshifts are 0.61 and 0.39. According 
to our previous discussion, the relative variability decreases 
with increasing luminosity, and, given the differences of the 
median values of the luminosities of the objects in these two 
classes, the real difference in the corrected relative variabil- 
ity of RL and RQ should be even larger. Actually, RL tend to 
be more variable than RQ over all magnitudes and redshifts, 
as can be seen in figure 12, where we plot median variability 
values calculated in bins with equal number of objects. Our 
observations then confirm earlier results obtained by Pica & 
Smith (1983). 



6 SUMMARY AND CONCLUSIONS 

In this paper we have presented the first results of a photo- 
metric monitoring of quasars that have been conducted at 
CNPq/Laboratorio Nacional de Astrofisica, in Brazil, since 
1993. The observations discussed here were done in the V 
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band and span the period between March 1993 and July 
1996. This set of observations has good photometric accu- 
racy (~ 0.019 V mag), achieved through differential pho- 
tometry with CCD detectors, what allows the detection of 
low levels of variability. 

The relative variability observed in the V band, 5 b s , 
decreases with increasing luminosity and redshift. However, 
our observations in the V band sample different rest-frame 
frequencies for quasars at different redshifts, and it is con- 
venient to convert all observations to a single rest-frame fre- 
quency (e.g., that of the V band) in order to compare the 
variability properties of these objects. In computing 8 res t, 
the relative variability expected in the rest-frame V band of 
the quasars, we have taken into account a model for the in- 
crease of the variability with frequency, assuming that an in- 
crease in luminosity corresponds to a hardening of the spec- 
trum. This model turned out to be consistent with the nu- 
clear variability observed during multi-wavelength observa- 
tions of the Seyfert galaxy NGC 4151 (Edelson et al. 1996). 

We have noticed that correcting the relative variabil- 
ity by applying the model above decreases the scatter in 
the anti-correlation of 5 with z or L, increasing slightly the 
significance of the relations between these quantities. Such 
anti-correlations are in agreement with results of other stud- 
ies, like Hook et al. (1994) and Cristiani et al. (1996). Note 
that the strong correlation in the Mv x z diagram for the 
objects in our sample preclude us of disentangling the real 
dependences of the variability on these two parameters. 

We have verified that the relative variability 8 re3 t tends 
to decrease with the luminosity in agreement with the ex- 
pectations of a simple Poissonian model, where the light 
curves are assumed to be the result of a random superpo- 
sition of pulses of same energy. Characteristic energies for 
individual pulses derived with such a model are consistent 
with those predicted by the nuclear starburst scenario (e.g., 
Cid Fernandes, Aretxaga & Terlevich 1996). Assuming that 
there are no differences in n, f var and t p for both classes, the 
pulse energies for the radio-loud objects are larger than for 
radio-quiet objects by almost one order of magnitude. The 
radio-loud objects in our sample tend to be more variable 
than radio-quiet objects. 
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